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ABSTRACT 

This paper shows that the plane wave expansion can be a useful tool in obtaining 
analytical solutions to infinite integrals over spherical Bessel functions and the deriva- 
tion of identites for these functions. The integrals are often used in nuclear scattering 
calculations, where an analytical result can provide an insight into the reaction mecha- 
nism. A technique is developed whereby an integral over several special functions which 
cannot be found in any standard integral table can be broken down into integrals that 
have existing analytical solutions. 
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1. Introduction 

The plane wave expansion, or the Rayleigh equation, is commonly used in nuclear 
physics; in particular it can be found in any standard text on scattering theory [1]. In 
the plane wave appproximation, a particle incident on a nucleus has a wavefunction 
approximated by a plane wave of the form e*^"^, where k is the particle's wavenumber 
and r is its position. This approximation is often used to get an insight into the reaction 
mechanism of the process studied, without the complications arising from using the 
more accurate distorted waves. 

Using the plane wave approximation will result in radial integrals that involve 
spherical Bessel functions. On the other hand, there are many other applications in 
nuclear physics that encounter these types of integrals amongst which are distorted 
wave calculations [2] and nuclear response function calculations [3]. These types of 
integrals have in the past been calculated analytically (many references exist, see for 
example references [2] and [4-11]), or numerically using complex- plane methods (see 
for example references [12-14]). In this paper, an additional advantage to the plane 
wave expansion is presented, whereby the expansion itself is used to derive analytical 
solutions to integrals over spherical Bessel functions and identities involving them. The 
technique itself is not new [11], but the author is not aware of any prior paper that has 
highlighted its versatility. In section 2, the plane wave expansion is introduced from 
which the integral representation of the spherical Bessel function is derived. Section 3 
builds on that by deriving integrals involving one spherical Bessel function. Section 4 
details integrals involving two spherical Bessel functions, amongst which is the Closure 
Relation for spherical Bessel functions. Section 5 deals with integrals involving three 
spherical Bessel functions. Section 6 extends the previous results to a larger number 
of any special or elementary function. Here the Closure relation is used to simplify an 
integral that cannot be found in any standard integral tables [15]. This can be done, 
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for example, by breaking an integral over 4 special functions into an integral over two 
integrals of 3 special functions each. Concluding remarks then follow. 



2. The Plane Wave Expansion 



The plane wave expansion, also known as the Rayleigh equation, is given by 



00 L 




4^E E i''yi'{r)Yf*{k)jUkr) 



(2.1) 



L=0 m=-L 



where y^(f) is the spherical harmonic function for the unit vector f, and jiikr) is the 
spherical Bessel function for A; > 0, which is assumed for the rest of the paper along 
with kn > 0, for integer n. An alternative formula is 



using the Vector Addition Theorem given by 

M=L 

= 2rTT ^ Yt'{r)Yt\k), (2.3) 

M=-L 

where Pl(cos^) is the Legendre polynomial and 6 is the angle between k and r. 

To limit the sum over L in eq. (2.1) or eq. (2.2) to a specific value, one needs to 
use the orthognality of the spherical harmonics 



oo 




(2.2) 




(2.4) 



where Q,t is the solid angle subtended by k. Alternatively, we can use the orthognality 
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of the Legendre polynomials 

1 

J Pl(cos d) Pl'{cos 6) d(cos 6) = ^^-^ 5l,l' (2.5) 
-1 

to arrive at the integral representation of the spherical Bessel function 



jL{kr) = j PLicosd)e"'''°'Uicosd). (2.6) 

-1 

Setting k = 0, one arrives at the identity 

jL{0) = 5l,o (2.7) 

3. Infinite Integrals Over One Sptierical Bessel Function 
Integrating eq. (2.6) over r and interchanging the order of integration we obtain 

J jL{kr)dr = j Pl(cos^) j e*^^^'^'^ drd(cos^). (3.1) 



-1 



However, 



oo 



j e'^^'dr = 2TTS{k), (3.2) 

— oo 

resulting in 

J 3L{kr)dT = i^7rPi(0). (3.3) 

— oo 

Another technique is to integrate the plane wave, eq. (2.1) over all space and using 

d\(^k.r ^ ^27rf5^{k), (3.4) 



one reaches 

oo 

J r'^jo{kr)dr= 2Tr'^S^{k), (3.5) 



which shows the spherical symmetry of S^{k), since there is no angular dependence on 
the left hand side. Also, since 

/pik . f 9„2 
(3.6) 

we can reach, using eqs. (2.1) and (2.4), 



TT 

jo{kr)dr = (3.7) 



4. Infinite Integrals Over Two Spherical Bessel Functions 
Again the plane wave expansion gives 

^i{k+k,).r ^ (4^)2 ^ ^ (i)^i+^^y^^i(^i)y^f^^(^2) 

Li, Ml L2,M2 (4.1) 

X yt:''ir)Y^2\r)jL,{kir)jLAk2r). 
Integrating over all space results in fixing L2 = Li = L and M2 = Mi = M 

{27Tf 6\ki + = (47r)2 ^(-l)^Fi^(fci)y/^*(fc2) 

L,M 

00 (4.2) 
X / jL{kir) jL{k2r)dr. 


Using 

5\k + h) = ^^^^^ S{n_^,^ -n^j, (4.3) 

multiplying both sides of eq. (4.2) by y^*(^i) ¥^(^2)5 followed by integrating over 



solid angles ki and /c2 and using 



Yi'i-k) = {-ifYifik), (4.4) 



results in 

oo 



j 3L{kir)jL{k2r)dr = ^^{ki - k2), (4.5) 
^ 

which is known as the Closure Relation for spherical Bessel Functions. 
Another relation that can be derived using eqs. (2.6) and (3.2) is 

r r/ 1 

r (A\L+L' r 

J jL{x)jL'{x)dx ==^—^ J PL{y)PL'{-y)dy. (4.6) 

— oo —1 

Now 

PLi-y) = (-l)^PL(y), (4.7) 

just as 

hi-x) = (-I)^Jl(x), (4.8) 

which reduces eq. (4.6) to 

oo 

J jL{x)3L'{x)dx = 2L^i ^L,L', (4.9) 

— oo 

known as the Orthogonality Relation for the spherical Bessel functions. Using this 
equation together with the plane wave expansion, eq. (2.2) leads to 



Pl{x)=^-^ I e'^'^jMdy, (4.10) 

— oo 

which is an integral representation for the Legendre polynomial valid for — 1 < a; < 1. 
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Also, dividing eq. (4.1) by and integrating over all space results in 

27r2 



|A;i + A;2| 



L,M 

00 

X j jLihr) jLihr) dr. 




(4.11) 



Now 



where A;< and /c> arc the smaller and larger of ki and k2, respectively. So substituting 
eq. (4.12) into eq. (4.11), multiplying by y^*(A;i) y^(fc2) and integrating over solid 
angles ki and k2 results in 

jL{hr)jUk2r)dr = 2(2ZTl) (f)^' ^^"^^^ 



which agrees with reference [4], equation 1, page 405. 

An identity can also be derived by equating the expansion in eq. (4.1) to the 
alternative expansion in terms of jl(|A;i + k2\r) resulting in 



(ifYi^ik,+k2)ju\h+k2\r) = v^Y.Y. ^^"'^'' r'^^'^T^i)^'^ 

Li,Mi L2,M2 V \ ^ ) 



X < L1L2OOILO >< L1L2M1 -M2\LM > Y£^^'{ki)Y-^^{k2) 
X jLiikir) jL2{k2r), 



(4.14) 



using 



'(2/1 + 1) (2/2 + 1) 



47r(2/ + l) (4.15) 

X < /1/2OOI/O >< /l Z2 "7-1 777-2 I / m >, 



where < li I2 mi 7x12 \ lm > is a Clebsch-Gordan coefficient, coupling angular momenta 
li and I2 to angular momentum I such that mi + m2 = m. 

If ^1 and k2 point in the same direction, i.e ki = k2 = k, then multiplying eq. 
(4.14) by Y^*(k), integrating over solid angle k then summing over Mi and M2 using 
the unitarity of the Clebsch-Gordan coefficients 

< Li L2 Ml M2 I L M >< Li L2 Ml M2 I L' M' > = 5l,l' Sm,M'^ (4-16) 

Mi,M2 

equation (4.14) reduces to the identity 

Mh + k2)r] = ^ < X1L2OOILO >2 

L,M (^^+i) (417) 

X jLiikir) jL2{k2r). 

Similarly, 

M{ki-k2)r] = (^)^-^-^ ^'^^Xyi'f.^^'^ < LiL200\L0 

L,M l^i^+lj (418) 

X jLiihr)jL2{k2r). 

Redefining x = kir and y = fc2r, multiplying eq. (4.18) by ji{x) and integrating over x 
would result in an equation of the form 

00 

J 3L{x-y)3i{x)dx = TT 5^ (i)'-^^^^ < L/OOIL2O >2 JlM, (4-19) 

-00 -^2 

using eq. (4.9) and the symmetry properties of the Clebsch-Gordan coefficients. 
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5. Integrals Over Three Spherical Bessel Functions 

Using the same procedure as before, integrating over the expansion of e*('^i+^2+fc3)-r 
and using eq. (4.15) amounts to 



s^{h+k2+ks) = 4 E E E (^)^' 



+1/2+^3 

Li,Mi L2,M2 L3,M3 



X 

oo 

X 





^ ^^'(2L3 + 1)^ < LiL200\L30> K{Li,L2,Ls,ki,k2,k3) (5.1) 



oo 



where 



K{Li,L2,L3,ki,k2,k3) = E E E < ^i^sMiMalLsMg > 

Afi Af2 M3 



(5.2) 



Multiplying both sides of eq. (5.1) by K*{L'^, L2, L'^, ki, k2, k3) and using the orthog- 
nality of the spherical harmonics 



(5.3) 



J dnj^^dnj^^dflf^^K*{L[, L'2, L'g, ki, k2, k) K{Li, L2, L3, k, h) 
= (2L3 + 1) 5l[,l, ^l'^m ^l'^M^ 

results in 



00 

<LiL2 00|L0> j r'^ jLiikir) jL2{k2r) jLaiksr) dr 


~ ~^ V(2^i + 1)(2^2 + 2)(2L3+Ty 

X j dVti^dVLi^dVt-^^ K*{Li, L2, L3, ki, k2, h) S^{ki + k2 + h)- 



Now 

5'{h + k2 + h) = ^(^3-|fel + fc2|) _ ^ ^_ 

Using 

,(,3-|M + fel) = ^*(cos«-^i±M), (5.6) 

where 6 is the angle between and k2, and the completeness relation for the Legendre 
polynomials 

or I 1 

^(^ - 2/) = E Pl{x) PM, (5.7) 



gives 



(53(^1 + ^2 + ^3) = -n.^)V(2L+l)PL(A)PL(cos^), (5.8) 



where 

A = 

2A;iA;2 



^2 I ^.2 _ 1,2 



and the factor /3(A) is given by [11] 

/5(A) = ^, A = ±1 

= 1, -1< A < 1 (5.10) 

= 0, otherwise. 

Equation (5.4) then becomes 

oo 

<LiL2 00|L0> j r'^ jLiihr) jL2{k2r) jLsihr) dr 


8hk2k3 y(2Li + l)(2L2 + 2)(2L3 + l) 

X f dl)j^^dl7^^dK*(Li,L2,L3,^i,fc2,-(^i + ^2))E(2^ + l)^i(^)^i(co«(^))- 

(5.11) 
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The integral in eq. (5.11) can be easily evaluated using eq. (2.3), eq. (4.4) and the 
Solid-Harmonic Addition Theorem for spherical harmonics [16] 

^3 J± I JZ /o r „ _i_ 1 \ 1/2 



L/ — U 

X (|)^E <iL3-L)LiMs-M)M\LsMs>Y^^rL''ik)Yi'ik2), 



M 

where the binomial coefficient is defined by 

jnj [n — m)lml 

The result of the integration is then the formula for the integral over three spherical 
Bessel functions 



< 



L1L2OOIL3O > J r'^ jLiihr) jL^ihr) jLaihr) dr 




4/C1/C2/C3 



X 



V < Li(L3-L)00|/0 >< L2LOOI/O > I i P;(A), 



Jl J2 J3 

where ^ ^ is a 6-j symbol defined in any standard angular momentum text 

[17]. An interesting special case of the above equation is when L3 = 0, which imposes 

Li = L2 = A leading to 

00 



and the inverse equation 

fei+fc2 2 2 2 

Mhr)Mhr) = ^^ I .,i„(MA(^^)*3. (5.16) 

|fc2-fcl| 
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6. Generalisation 



It is possible to extend the results of the last section to finding analytical solutions 
to infinite integrals over several special or elementary functions (provided the integral 
exists) if the integral over a fewer number of these functions combined with a spherical 
Bessel function is known. One can show, using the Closure Relation for spherical Bessel 
functions eq. (4.5), that 



oo oo 

j r'^ XLi{ki,r) XL2{k2,r) XL3{k3,r) XL4{k4,r)dr = ^ J k^dk 



(1 \ (1 

X I / r'^XLAki,r)xL2{.k2,r)jL{kr)dr\ I r''^ XLzik^^r') xiAk^^r') jL{kr')d\ 



'^V I , 



(6.1) 

where XL{k,r) is a special function of order L that, in general, depends on k and r. 
So, in effect, the integral over four special has been broken up into two integrals over 
3 special functions. The final integral over k may be easier to handle than the initial 
integral, or in the least is an alternative approach to solving the integral analytically. 
The following example illustrates the advantage: 

Example- Evaluate the following integral over 4 spherical Bessel functions analyti- 
cally 

oo 

^ = J f'^hikr) ji{kr) j2{kr) j2{kr)dr. (6.2) 



Solution- The integral can be rewritten as 

oo / oo \ / oo \ 

I = - q^dq / r'^ji{kr) ji{kr) jo{qr)dr / r''^j2{kr')j2{kr')jo{qr')dr' . 

(6.3) 
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Using eq. (5.16) this reduces to 



2k 


This can be easily evaluated using Pi{x) = x and P2{x) = ^ (3a;2 — 1) to give 
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7. Conclusions 



The plane wave expansion was utilised to derive integrals over one, two, and three 
spherical Bessel functions. The generalisation to a larger number of special functions is 
exhibited through the use of the Closure Relation for spherical Bessel functions, which 
is also derived here using the plane wave expansion. In addition, several identities 
involving spherical Bessel functions were also derived using the same techniques. 
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